Incipience of quantum chaos in the Jahn- Teller model 
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We studied complex spectra of a two-level electron system coupled to two phonon (vibron) modes 
represented by the E®e Jahn- Teller model. For particular rotation quantum numbers we found a 
coexistence of up to three regions of the spectra, (i) the dimerized region of long-range ordered (ex- 
tended) pairs of oscillating levels, (ii) the short-range-ordered (localized) "kink lattice" of avoiding 
levels, and (iii) the intermediate region of kink nucleation with variable range of ordering. This 
structure appears above a certain critical line as a function of interaction strength. The level clus- 
tering and level avoiding generic patterns reflect themselves in several intermittent regions between 
up to three branches of spectral entropies. Linear scaling behavior of the widths of level curvature 
probability distributions provides the conventionally adopted indication for the presence of quantum 
chaos. Level spacing probability distributions show peculiarities of the partial (for fixed quantum 
angular momentum) as well as of the cumulative (all angular momenta) case. The clustering of 
levels with two and three dominant spacings at fixed angular momenta causes notable deviations of 
the cumulative distributions from the Poissonian one. 

PACS numbers: 05.45.Mt,31.30.-i,63.22.+m 



I. INTRODUCTION 

Generally, a multitude of avoided crossings in complex 
spectra of a quantum system is considered as the signa- 
ture of chaos which appears in a quasiclassical limit of 
the respective nonlinear many-body system with repul- 
sive interactions [3, 0. In solid-state physics, quantum 
chemistry, or quantum optics investigations of complex 
excited energy spectra of models with two electron lev- 
els coupled with one or more boson modes are of special 
relevance for understanding experimental - e.g., optical 
or transport-properties. 

Numerical investigations of the excited spectra of cer- 
tain two-level one- and two-boson systems fe.g., exciton 
@, Q and E®e Jahn-Teller (JT) model [|] show the 
existence of level avoidings and related highly excited 
"exotic" (localized) states. Signatures of chaos in two- 
level boson systems remain objects of interest since the 
early studies by Hamiltonian methods [f| 0, H, S EH El • 
Later on statistical methods based on the random matrix 
theory - e.g., level spacing and level curvature probabil- 
ity distributions - were applied mostly for one-phonon- 
mode models. The level spacing distributions show ir- 
regular behavior between the limiting Poisson (regular) 
and Wigner (fully chaotic) distributions [HI, EH, 0, EH • 
Specifically, among the level spa cing distributions of the 
two-level one-phonon model [13, EH Ej] there appeared 
M-shaped distributions with two symmetric peaks for two 
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dominant spacings. Such distributions are a signature of 
level clustering (dimerization). 

The most familiar representative of the two-level two- 
phonon models is the E<§5e JT model [l6| with linear 
coupling to phonon (vibron) modes of different parity 
against reflection. This model is a prototype for phonons 
removing the degeneracy of electron levels by an inter- 
action with the antisymmetric phonon mode and tunnel- 
ing between the levels assisted by the symmetric phonon 
mode (in the one-phonon case the tunneling is purely res- 
onant). The presence of two-phonon modes with equal 
coupling strengths imposes additional rotational symme- 
try to the model. Differences between the one-mode 
two-level phonon model and two-mode JT models be- 
come evident when one compares their dimensionality 
after elimination of the electron (level) degrees of free- 
dom: while in the first case the reduction results in a 
one-dimensional nonlinear (with self-interaction) quan- 
tum oscillator, in the two-mode JT cases the elimination 
yields a two-dimensional quantum oscillator with nonlin- 
ear coupling of its components. While in the one-boson- 
mode case one finds correlations of two sets of excited 
states, in the two-boson-mode case mixing of up to three 
sets takes place depending on the range of parameters. 
As we show, this leads to the appearance of a third peak 
in the level spacing distribution (trimerization) for the 
respective range of parameters. 

Recently, Yamasaki et al. [17] for the first time inves- 
tigated the possibility of chaos in spectra of the E(g>e 
JT model. Their analysis was based on the approxi- 
mation of the Hamiltonian by the adiabatic "Mexican 
hat" potential. Additionally the Hamiltonian was sup- 
plemented by an explicit term with nonlinear mode cou- 
pling of a trigonal symmetry in order to simulate the 
effects of fluctuations and nonintegrability. Thus, non- 
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linearity was included via mode coupling in addition to 
the mean-field bare part of the Hamiltonian. The au- 
thors concluded that quantum chaos reflects itself in the 
Wigner-type level spacing distribution as a consequence 
of the said nonlinearity of the Hamiltonian. Meanwhile, 
for the linear part an absence of such patterns was stated. 
Those conclusions were confirmed at the classical level as 
well. It is to be noted here that the way of passing to 
a semiclassical approximation in spin-boson systems is 
not unique and presents an essential ambiguity from dif- 
ferent possible ways of decoupling [6| . Different ways of 
performing a semiclassical approximation are known to 
lead to different answers concerning the chaotic behav- 
ior of the system. This ambiguity means that a classical 
analog of such a model is not well defined and cannot be 
a reliable object in exploring the quantum chaos issues. 

In this paper we investigate thecharacteristics of ex- 
cited (quasi continuum) spectra of E®e JT model. The 
present approach differs from that by Yamasaki et al. in 
the following: (a) we do not introduce explicit nonlin- 
earity into the initial Hamiltonian. However, its SU(2) 
symmetry involves an intrinsic nonlinearity which is re- 
vealed by exact elimination of the electronic degrees of 
freedom [l8|. (b) We start from the E<g>(6i + b 2 ) model 
with different coupling strengths for both modes. The 
rotation-symmetric E®e model represents its particular 
case with equal coupling constants and thus with the 
symmetry of a higher (rotational) symmetry group (the 
difference of the coupling constants in realistic systems 
is likely to be caused, for example, by spatial anisotropy 
of crystals). The importance of starting from a more 
general situation is evident from the nature of quan- 
tum fluctuations. Namely, the variational approach to 
the E(g>(&i + £2) model used for the calculation of the 
ground state [18| yielded the largest deviations just for 
the rotation-symmetric case. In the ground state the 
abrupt change of energy at equal couplings is an artefact 
of the adiabatic approximation: the energy region of the 
quantum E(g>e model is situated within a smeared border 
around (3/ a = 1 (a and (3 being the coupling strengths of 
the antisymmetric and symmetric modes, respectively) 
between the " selftrapping" (a > (3) and "tunneling" 
(a < (3) part of the ground state of the model with bro- 
ken rotational symmetry. In the smeared transition re- 
gion both phases coexist and are quantum correlated (en- 
tangled) via nonlinear correlations between the modes. 
There the phonon-assisted (by the symmetric mode) tun- 
neling contribution to the energy (from the admixture of 
two reflection symmetric levels in the excitation reflec- 
tion Ansatz for the variational wave function [18|) causes 
an essential decrease of the ground state energy because 
of changing parity of the wave function by the reflection 
operator and results in the vanishing of the Ham factor. 
The phase transition at a = f3 to the rotational symmetry 
phase provides a representation by the additional rota- 
tional quantum number j. We will show that at a given 
j in Eg)e model analogous phonon assisted tunneling be- 
tween adjacent excited levels occurs and is responsible for 




FIG. 1: Energy spectra for j — 19/2,255/2 as a function of 
a. 



the flip (kink) between the levels [this mechanism can be 
well recognized, e.g., in Fig. 1(b)]. This nonlinear effect 
causes a mixture of the adjacent (even and odd) levels 
at stronger couplings and leads to nucleation of the new 
kink lattice phase when increasing the coupling constant. 
Naturally, the tunneling mechanism due to short-range 
quantum fluctuations at large j is then the origin of the 
peculiar irregular behavior of related level spacing distri- 
butions and does not allow for full developing of Wigner 
chaos. The phonon-assisted tunneling and nucleation of 
the kink lattice phase can also be demonstrated within 
the formalism of a Calogero-Moser gas of pseudoparticles 
with repulsive interaction [l9|, [2(| ■ 

In Sec. II we briefly present an analytical formulation 
of the eigenvalue problem in radial coordinates. Since 
the literature on one-phonon models is more abundant 
and certain similarities between both models appear, it 
is useful to point out the essential differences based on a 
comparison of the relevant integrable models. In the rep- 
resentative one-boson mode model considered by Cibils 
et al. 13] the integrable case was labeled by a pair of 
quantum numbers //, v where /i = — s, — s + 1, . . . , +s and 
v = 0, 1, 2 In our model the integrable case is labeled 
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by three quantum numbers n r ,j,p, where n r is the ra- 
dial quantum number n r = 0,1,2,..., j is the rotation 
quantum number j = 1/2,3/2,5/2,... and p is parity, 
p = ±1 for gerade and ungerade states. Thus, there are 
two groups of states which represent a doubly degenerate 
set characterized by parity p = ±1. Evidently, n r is anal- 
ogous to v, p to s = ±1/2, and j has no analog because 
our oscillator is two dimensional. The Hamiltonian anal- 
ysis is supported by numerical evaluations of the excited 
energy spectra and corresponding phonon wave functions 
in a representation with definite parity and square of 
angular momentum. We choose a representation of the 
spectra and their probabilistic characteristics (Sec. IV) 
for definite j because it contains well-defined physical in- 
formation. However, the real spectra and probabilistic 
distributions are composed of the contributions of all j's. 

In Sec. Ill we present approximate analytical ap- 
proaches to the E<g>e JT model for the strong and weak 
intcrcluster coupling in order to elucidate the peculiar 
features of the numerical excited spectra and wave func- 
tions. In particular, we find the approximate analyti- 
cal form of wave functions for strongly localized ("ex- 
otic") states as well as for extended states and compare 
them with respective numerical results. Special interest 
is given to identification of the regions of the spectra with 
different extent of ordering, from the long-range ordered 
region of linear oscillations (extended states) to the short- 
range ordered region of localized states with dynamical 
avoiding (kink lattice). The dynamical Calogero-Moser 
approach appears to be capable of describing all three re- 
gions of the spectra: the dimerization region with pairs 
of oscillating levels, the kink train lattice, and the inter- 
mediate region of nonlinear fluctuations accounting for 
the kink nucleation phase. The intermediate region was 
found to exhibit quantum chaotic behavior under cer- 
tain conditions. Section IV complements the dynami- 
cal insight by that based on statistical properties of the 
spectra. We consider the nearest level-spacing proba- 
bility distributions and probability distributions of level 
curvatures used as relevant statistical characteristics of 
complex spectra [l|, d, [2l[ . The particular level-spacing 
distributions for fixed rotation quantum numbers are be- 
lieved to be the most informative, since they behave very 
specifically for different parts of the spectra, reflecting 
their specific dynamic behavior. 



II. HAMILTONIAN AND EQUATIONS FOR 
EXCITED SPECTRUM 



We investigate local spinless double degenerate elec- 
tron states linearly coupled to two intramolecular phonon 
(vibron) modes described by Hamiltonian 

H = (b\b 1 +blb 2 + l)I + a(b\+b 1 )a z -(3(bl + b 2 )a x , (1) 



where a T — 



0, 1 



i. i)/' ~ \ i,o y ° z ~ 1^0,-1 

are Pauli matrices, I is the unit matrix. This pseudospin 



0,-1 



notation refers to two-level one-electron (spinless) sys- 
tem. 

For general a ^ (3 the Hamiltonian possesses a special 
SU(2) symmetry with the reflection operator R acting 
on phonon and electron subspace: 



R = R„iR 



■ph • 



R, 



cr x , R P h = exp(i7r6|6i) (2) 



and RphQi = -QiRph, R P hQ2 = QiR v h\ i-e., phonon 
1 is antisymmetric and phonon 2 symmetric against the 
reflection. Then, [R, H] = and the eigenstates can 
be cast as eigenfunctions of the reflection operator with 
eigenvalues ±1 which are good quantum numbers. 

If a = f3, the Hamiltonian becomes rotationally sym- 
metric in the plane (Qi x Q 2 ) and there arises another 
good quantum number, the eigenvalue of angular mo- 
mentum J which in fact is nothing but the infinitesimal 
generator of the rotational group: 



J = - b+b 2 ) 



2<V 



(3) 



Therefore, in the vicinity of the symmetric JT case it 
is feasible to cast the solutions in the form respecting 
this rotational symmetry. To account accurately for the 
group properties of the Hamiltonian one introduces radial 
coordinates in the Qi x Q 2 plane Qi = rcos<p, Q 2 = 
r sin 4>. The angular momentum operator §3§ is 

* . d 1 

and the phonon part R p h of the reflection operator © 
acts as Rph(r,(j))f(r,<f>) = f(r, n - cj>) on some f(r,(f>) 
affecting only the <j> coordinate. The eigenfunctions of 
the angular momentum operator with eigenvalue j (i.e., 
J\]/ = jif>) can be set as follows: 



/i(r) ■ |+)e^-4) + f 2 (r) ■ |-)e^( i+ *), (4) 



with arbitrary fx (r) , f 2 (r) to be determined below. Elec- 
tronic wavefunctions are expressed as 



2U-! 



If 1-i 
2 I 1 + i 



The angular quantum number j takes on the values j = 
±1/2, ±3/2, ±5/2, . . .. The operator J does not commute 
with the reflection operator R, but J 2 does, so it is more 
convenient to build solutions as joint eigenfunctions of R 
and J 2 rather than J. 

In order to make use of the reflection property of the 
Hamiltonian we perform the custom Fulton-Gouterman 
transformation [22| by means of the operator U: 



U= — ( 

y/2 \ l,-R P h 



(5) 
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The resulting Hamiltonian is diagonal in electronic space 
and is represented in its full form as follows: 



H = UHU- 1 = H r 



2^ ' d6 



V2ar ( C0S ^ _ S ^ n < ^^'P h ^ 

I cos 4> + sin 4>R p h 

+ (a-/3)V2r.f sin ^ . ° 

\ -sin 0i? p h 



where 



I f d\ 1 2 



(6) 



(7) 



is the radial part and R p h , Eq. ([2]) is the (highly nonlinear) 
phonon reflection operator in the radial coordinates. In 
Eq. (J6j> we extracted explicitly the last term ~ (a — (3) 
which breaks the rotational symmetry of the problem. 
Just for reference we mention the corresponding trans- 
formation of the angular momentum and its square: 



T d 1 o 

= ~ l d6 ax + 2 av ph 



J 2 = - 



9 r? 
■jr-cr z K p h- 



(8) 



(9) 



The transformed electron reflection operator is of course, 
trivial per constructionem, R e i — a z . 

Although the E(8e JT model (a = (3) was treated for 
the first time by Longuet-Higgins et al. |23j a long time 
ago, we present briefly the solution using another repre- 
sentation of the states. This representation enables us to 
distinguish the extended and "exotic" (localized) states, 
as we explain in Sec. III. The case a ^ f3 can be included 
on equal footing. The wave functions we are concerning 
with have definite parity ±1 and square of angular mo- 
mentum j 2 . Let us denote K — exp — 1/2)) = ±1. 

i) In the case K = +1 (j - 1/2 = 
0, 2, 4, 6, . . . ; — 2, — 4, . . .) the gerade solutions read 



* = 



E 



Mr) 
2tt 



sin j 



1 



(10) 



cos( j + -(-!)' 



We shall refer to the angular wave functions in (jit 
|±> = ^[sin ((j + ±(-1)00) T cos ((j + U-lfm- 

ii) In the case K = —1 (j — 1/2 = 
1, 3, 5, 7, ... ; — 1, —3, . . .) we have the expression similar 
to (JTUJ) but with the + sign in the square bracket. 

The ungerade case can be recovered if one takes the 

'oS . , , (i 



vector 



1 



instead of 



We limit ourselves to the 



case of gerade parity that is, in the representation ([6]) we 
are to take the upper row of the matrices. 

The secular equations for the vibronic part (functions 
/i( r ),/2(f)) in both cases are 



H 'fi+2?2 (i+o(-i)" 



fi+aV2rf k = Efi, i ^ k, i,k 

(11) 

where the radial part H r is given by frj). 

In the case a — (3 — the problem is exactly that of 
two independent oscillators and its solution in terms of 
special functions is well known. Thus, the solutions are 
grouped according to the |j| value (j — 1/2, 3/2, 5/2, . . .; 
the parity is +1, as before). Within the same j the func- 
tions are divided into two groups, 

(a) States with f% = Q and 



= 1,2 



* + K, j;r) 



C 



:exp (_L.). r li-l/2| . lFl 



-n r , 1 + 



J -2 



where n r ("radial" quantum number) takes on positive 
integer values n r = 0, 1, 2, . . .; i-Fi(a, b, z) is a confluent 
hypergeometric function [satisfying equation zy" + (b — 
z)y' — ay = 0], and the norm is 



1 



n, + b-l/2| 



(13) 



li-i/2| 



(here = n\/[m\(n — m)\] is a common combina- 
torial factor). We label these states by + as they 
correspond to the |+) -function of the angular momen- 
tum. The product of two radial functions is defined 

oo 

as J r<f>(n r , j, r)$(n r , j, r)dr. The energy spectrum is 
o 

E + = 2n r + \j — 1/2 1 + 1, and the "main" quantum 
number is therefore n = 2n r + \ j — 1/2|. 

(b) States with /i=0 and f% = $(~>(n r ,j;r) and en- 
ergy E~ = 2n r + |i + 5 1 + 1 are given by the same ex- 
pression as (TT21 but with j + 1/2 instead of j — 1/2. We 
refer to these states as |— ) states. 

The two groups mentioned correspond to two initially 
degenerate electron levels; these states are coupled to- 
gether by means of the strongly nonlinear transforma- 
tion involving the reflection operator R(r,4>)', thus, the 
situation is similar to the Foulton-Gouterman treatment 
in ordinary phonon coordinates, where the nonlinear re- 
flection operator coupled the phonon states pertaining to 
the lower and upper electron level [HI]. 

Writing secular equations for the problem with a =/= 
is straightforward. The solution is sought as f± = 
J2n r a n r ®i r > h = Et, r & Jir $ n r ) or > in unified notation, 

oo 

f = H c„<&„, where {c„} = {a , bo, ai, b\, a 2 , b 2 , a 3 , . . .}, 

n=0 



(12) 
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and 



$ — 



$ + (n r , j, ?*), n = 2n r , 
$-(n r , j,r), n = 2n r + l. 



(14) 



The only non-zero off-diagonal elements arc 



f nrnr = V^J *+ $" r r • rdr = n r + 1 + j - ^ (15) 

/„ r+ i„ r = y/2 J $+ r ■ rdr = -s/2-J^~F{lQ) 

Therefore, the secular equation acquires the familiar 
tridiagonal form (E^ is "unperturbed" energy for the 
corresponding n,j), 

C n E^ + a(c n+ if nrnr + Cn-lfnrtir-l) = Ec n, " = 2n, 
C n E^ + a(c n -lfn r n r + C n+ lfn r + ln r ) = Ec n , 71 = 2n r + 

with off-diagonal coefficients given by Eqs. (fT5)) and (fT())) . 
The asymmetry of the coefficients with increasing j is 
crucial for the appearance of nonlinear effects in spectral 
properties of the model. We remark by passing that Eqs. 
(fl~5|) and (fl6|) differ from those given in [23| by the "mi- 
nus" sign near f n +in- Changing this sign to plus means 
merely redefinition of the base wave functions (every even 
base function acquires the opposite sign) and does not af- 
fect the energy spectrum. In the following it will be more 
convenient to use the + sign of both expressions (fl~5|) and 

(USD. 

By numerical diagonalization of the system (|17p and 
(|18[) we obtain the phonon excited energy spectrum as 
function of a (we have taken n up to ~ 1200), as well as 
the corresponding wave functions. Figure 1 shows exam- 
ples of the spectrum as function of a for j = 19/2, 255/2). 
The sample numerical wave functions x-n of a few excited 
states of (P) are shown in Figs. 2 and 3. Figure 2 presents 
the wave functions \n in phonon coordinate representa- 
tion (in the plane Q\ x Q 2 ), while Fig. 3 shows the pro- 
jection Xn, m = KXnl^ro)! 2 °f the exac t numerical wave 
functions on the unperturbed (a = 0) states of the radial 
oscillator (JT^J) [Hj]. The right-hand side of Fig. 4 visual- 
izes the whole picture of wavefunctions at a given a, j by 
assigning to each state C^ n ' = {Ci, C2, . . . , C m , . . .} the 
spectral entropy S — — J2 m l°g Cm (which gives the 
width of the wavefunction in the representation of Fig. 

One remarkable feature of excited states in JT mod- 
els (both rotation symmetric and that with broken rota- 
tion symmetry) are strongly localized ("exotic" [3j) states 
emerging abruptly between more extended ones. These 
states are exemplified in Fig. 2 (n = 61,63) and in Fig. 
3 (n = 58 for j = 19/2). The localization of those states 
is remarkable in comparison with neighboring or even 
lower-lying ones (like n = 38, 62 from Fig. 2 compared to 
n = 61,63). 

There is an intimate relation between the position of 
localized states and spectral picture of levels (Fig.l). 



Closer inspection [see the insets of Figs. 1(a), 1(b)] shows 
that each and every energy level taken for a given j is 
avoided. Varying the parameter a infinitesimally slowly 
will cause the system to follow one and the same adia- 
batic level. But on larger scale the spectra of Fig.l show 
two sets of "diabatic" lines 0] (directed upwards and 
downwards) which reflect a different behavior of the two 
main wells of the effective potential [18j when changing 
a. From the right-handside of Fig. 4 it is possible to note 
that the states for large a, j are grouped into two sepa- 
rated branches and the domain of localized states with 
small spectral entropy emerges at a critical energy cor- 
responding to the appearance of the upper sheet of the 
effective potential. If one varies the parameter a fast 
enough the system will follow the " diabatic" lines chang- 
ing the energy level number so that its wave function 
^^mains almost untouched which means that the system 
resides on the corresponding sheet of effective potential. 

"exotic" states appear to pertain to the directed up- 
wards diabatic lines of the spectrum seen markedly in 
Fig. 1. The multitude of avoided crossings (and there- 
fore the domain referred to as " quantum chaotic" ) is thus 
formed in the intermediate energy region where these two 
"phases" coexist- that is above the first (lowest) diabatic 
line. It is clearly seen in Fig. 1(b) for large values of the 
angular momentum number j, while for small values of 
j these exotic states, although still existing, do not come 
into complex interplay with the ordered states. This pic- 
ture is in agreement with the semiclassical treatment of 
the problem where the classically chaotic region is known 
to emerge sharply above some critical energy. 



III. ANALYTICAL QUANTUM TREATMENT 

A. Strong-coupling approximation: "Exotic" states 

The peculiarities of the spectra shown in Fig. 1 as well 
as their connection to the character of pertaining wave 
functions (Figs. 2 and 3) can be understood from merely 
a rough approximation to the set of equations (| 1 T[) and 
(fTSf in the following way. From Eqs. (| 1 5|) and (|16p . 
for \j\ > n r , there holds \fn r ,n r \ > \fn r +i,n r \> and the 
basis functions in Eqs. (fTT|) and (fT5j) appear to be "clus- 
tered" in the sense that the coefficients C2n r andC2n r +i 
pertaining to the same n r are coupled together (through 
fn r .n r ) essentially more strongly than to neighboring el- 
ements. [In order to avoid confusion we note that the 
terms "strong-" and "weak-coupling-approximation" for 
this section are not meant in the usual sense of large and 
small values of the parameter a but refer to the relative 
coupling inside and between the "clusters" (C2„, C-2n+\))- 
They are determined by the relation between f nn and 
fn+in~ i- e -, by the relative values of quantum numbers j 
and n]. Therefore, for the zeroth approximation of strong 
coupling f n+ i rl = 0, and the reduced version of Eqs. (fT7|) 
and (fT5|) for each n r gives us two sets of localized solu- 



FIG. 2: (Color online) Numerical wave functions \ n in the plane Qi x Q2, a = j3 — 2. States n=61 and 63 are localized 
("exotic"). 



tions for the energy 



(HHJ) then yields 



E nr j (a) = 2n r +\j-l/2\+3/2± v / 2(n r + 1 + | j - l/2\)a , 

(19) 

with corresponding localized eigenvectors C2n r +i — C 2rir 
(+ sign) and C 2 n r +i — -C 2rir (— sign) , the other C n 
being zero. From Fig. 1 we can conclude that the slopes 
of the corresponding lines follow the predictions given 
by Eq . (fl"9j) . This is valid for several lower clusters with 
small n r which show up as directed upwards and down- 
wards "diabatic" lines. At higher n r , when j < n r the 
intercluster f nr +in T and inside-cluster elements |/ ni .n r | 
become comparable and the long-range correlations yield 
the extended wavelike solutions seen in the upper parts 
of spectra in Figs. 1(a) and 1(b) [cf. inset of Fig. 1(a)]. 

To find the form of localized wavefunctions along the 
diabatic line we proceed by treating terms with f n +in as 
perturbation to the zeroth-order solution (fT!)]) . This ap- 
proximation is valid far from the points of avoided cross- 
ing where the wave functions are essentially localized and 
unaffected by neighbors. Let us denote SE = E — E* n 
where E^ is the zeroth-order strong-coupling solution 
(I19p for a given m and use the strong-coupling ansatz 
C = Ci, . . . , C 2n = C 2 „+i- The sum of Eqs.([T7D and 



n r = 0, 1, . . . 
[2(n-m)+a(f„ 



V2 



-/mm)]C , 2ll H — —a (Vn + lC 2 „+2 + VnC 2n -i) 

Further, let us formally mark the cluster (C*2n, C 2 n+i) 
as \n) and introduce formal creation and annihilation bo- 
son operators b + and b acting on these clusters. Since 
j ^> n, the term /„„ — f mm can be simplified by keeping 
only the leading-order term ~ [n— m). Then the effective 
Hamiltonian equation in strong-coupling approximation 
is (n|n) is replaced by b + b\n}) 



61 



ni? f b+b- 



-Lo ( 

V2 



{8E n + d£)m)\n) 



with 



n 



(+) 

eff 



2 + a/y/2j-l. 



(21) 



(22) 



The resulting system presents a quantum oscillator dis- 
placed by the term aj\[2. The effect of nonlinear fluc- 
tuations is reflected in the squeezing of its effective fre- 
quency (|22j) due to j down to its minimum f^// ~~ * ^ 
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n = 37,39,41,43,45,48,...; such exact wave function 
is exemplified in Fig. 3 for n = 58). We performed 
an extended comparison of the exact wave functions to 
those found via this method in a wide range of parame- 
ters a and j. Our calculations show that the solutions of 
Eq. (|2ip yield a very plausible approximation to the exact 
solutions of Eqs. (fTT)) and (TT8"]) for localized states, espe- 
cially for big j: for example, for several lowest diabatic 
states the projection (^l 1 !/) of the approximate over the 
exact state for j = 255/2 lies within 0.999 - 0.995. The 
same method is also suitable for approximate solutions 
found for several strong-coupled extended states, but 
there the cluster ansatz is antisymmetric, C^+i = — C2n- 
The corresponding equation differs from Eq. (|2"0|) by the 
— sign near both terms with a, respectively; there ap- 
pear the — sign of a in Eq. (|21[) and in the definition of 
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(23) 
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FIG. 4: Reduced energies (see text) and spectral entropies of 
wave functions Xn for different couplings a and j. 



in the quantum limit j — > oo. The term with m in fact 
labels the "sites" of a pseudolattice. In this approxima- 
tion the solutions for different m are translation invari- 
ant. The subsequent solutions of Eq. (|2"Tj) are obtained in 
a standard fashion applying the displacement operator 
D(i) = exp( 7 (fe+ - 6)), 7 = -a/(V^n[f f ) on the clus- 
tered pseudostates \n) = {C-z n , C^n+i)^ giving the local- 
ized functions on subsequent diabatic lines (for example, 
for j = 255/2 and a — 2 the lowest states of this ficti- 
tious oscillator correspond to wave functions with levels 



which is now antisqueezed due to nonlinear fluctuations. 
In the limit j — > oo both oscillators degenerate at fre- 
quency 2 and so the intercluster interaction effectively 
vanishes. In the absence of the intercluster interaction 
the energy lines from different clusters would regularly 
cross. The intercluster interaction leads to the spreading 
of the localized wavefunctions as described above and to 
the splitting of levels (avoided crossings). In the vicinity 
of an avoiding crossing we can approximately represent 
the relevant Hamiltonian as an effective (2 x 2) matrix: 



H = 



pr u 
u 



(24) 



where r = a — olq, is the point of level collision if 
there were no interaction. The matrix (|24p accounts for 
the collision of two levels, one lying at E = (if r = ±oo) 
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and the other moving asymptotically with velocity p with 
"pseudotime" r. In the representation (j2"4")l we are in the 
frame of reference moving downwards with a so that the 
velocity of the nonexotic levels is zero. If the colliding 
levels are weakly affected by neighbors we can consider 
p and u to be approximately constants. The minimal 
avoided crossing at r = is 2\u\. Including other non- 
exotic levels (thus situated at distances 1,2,.. .) leads to 
the extension of the matrix (|24[) to the shape given by 
Eq. (2) of the paper by Gaspard et al. [25| and account- 
ing for a single level moving with velocity p and crossing 
subsequently parallel energy levels y\ , y 2 , . . . with effec- 
tive interaction strengths u\, u 2 , ■ ■ ■■ It was shown there 
that for equidistant energies yi and equal strengths u% 
this system immediately gives a solution in the form of a 
soliton propagating with pseudotime r through the "lat- 
tice" y n with natural analytical predictions as to the pro- 
file of energy eigenvalues x n (r). This soliton like solu- 
tion corresponds to our diabatic level. The parameter 
p equals the relative slope of exotic (diabatic) and ordi- 
nary level and can be estimated in the strong-coupling 
approximation using Eq. (|19[) ; for the first diabatic line, 
it therefore holds p ~ 2y / 2j. In the same approxima- 
tion the level splitting (parameter u) can be estimated 
as (^ n \H(ao)\^ m ), where H is the exact Hamiltonian 
and ty n , $ m are wave functions in the strong-coupling 
approximation of the localized and extended states with 
numbers n and to, n = 1,2..., m = n + 1, . . . taken in 
the point ao of crossing of the corresponding unperturbed 
levels. For example, for the first avoided crossing of the 
lowest diabatic line the estimation gives u ~ 2j \fj. For 
large numbers to of subsequent collisions the gap value 
turns out to be approximately constant which suggests 
the applicability of the soliton like solution like those 
considered by Gaspard et al. (25j . 

B. Weak-coupling approximation 

This is the case opposite to the previous one and valid 
for j <C n r . Rewrite the system (fT?} and (| with coef- 
ficients (TT5l) and (fTT)|) with positive signs in the following 
way (here we again label by n instead of n r ): 



with form al notations: C 2n — > \2n) and C2 n +i — > |2n+l) 
and operators b + ,b acting upon our radial functions \n) 
as ordinary creation and annihilation operators [in the 
previous subsection the corresponding " pseudostates" re- 
ferred to the cluster (C 2n , C2n+l)]j the perturbation op- 
erator SV acts according to the recipe: 

5v{ |2n> -ok ./l 2n+1 > (28) 
° V \ \2n+l) ™ \ \2n) [Z *> 

The interaction (|28[) ties together the elements inside the 
cluster (6*2,1, CVn+i) and when it becomes considerable 
we again recover the strong-coupling approximation of 
the previous section. 

From Eq. ([2"7|) we immediately get the solution of the 
unperturbed problem Hoip = E^ip in terms of the 
displaced Fock states by acting the coherent operator 
D(^f) = exp(7(6+ — b)) on our pseudostates: 

|*) = \n) = D(-a)\n > (29) 

with unperturbed energies E~^ = n — a 2 . Therefore it 
is reasonable to investigate the reduced spectrum from 
which the "secular parts" of energies are subtracted: 
Sx n = E n — n + a 2 — (j — 1/2). On the left-hand side 
of Fig. 4 we plot the reduced wavelike spectrum 5 
function of n for typical values of a,j. Examining Fig. 4 
we observe that extracting the secular part reveals pat- 
terns which have not been noticed on the normal scale of 
the wavelike spectrum. The most interesting peculiarity 
is thus a crossover from two twisting quasi-sinusoids for 
large n and/or small a to the pictures which curiously 
resemble chaotic intermittent patterns containing rather 
characteristic windows of a regular motion. The cutoffs 
on the left parts of the spectra (small n) correspond to 
the first diabatic line below which such a reduction looses 
its sense since the perturbation SV can no more be con- 
sidered as small compared to Hq. 

The standard perturbation scheme applied to l[2"7j) 
gives up to the first perturbation order the expression 
for the spectral energies the notation h merely reminds 
one of the fact that the base functions are now "displaced 
Fock states" (f2"9")l not to be confused with \n)] 



, , . En =n - a 2 + SVnn + ■■■ (30) 

C 2n (E° n -E)+a [C 2n +i V^ T+T + C 2n -iV2 ^\ + aC 2n+1 k n = 0, 

C 2n +i{El n+1 -E)+a [C 2n V2lT+T + C 2n+2 V2lTT2] + S^kK = 



(25) 

with 



5Vnn = Y,{n\h){h\5V\s)(s\ 



h , s 



k n = V2^+l ( \/ 1+ 2^i _1 ) > n = M>--- (26) =2a^{2h+l\D{-a)\n)-{2h\D{-a)\n)-k h , (31) 



h=0 



(k n ~ j/v^n+T « 1 if j « n). The whole system ^ wher f K , ZS™? by , ® \ , ThC , elemcnts 

, - , (n to) = (n \D{6) to) can be calculated directly: 

is elegantly cast as iJ* = EV with \ \ / \ \ v/\ / J 

I m/mi \ ( P 2 \ a\n-m\ sgn^n-m) \ n - m \ , fl2 x 
(n\D(p)\m) = exp -— \fi 1 -==—L min{ ) {l3 ), 

H = b+b + a (b+ +b\+SV = H + SV (27) ( 32 ) 
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where L b a (x) are Laguerre polynomials. 

In the scope of the present paper we do not present 
further detailed analysis of Eq. ([3"Tj) . just mentioning now 
that it gives essentially the same structure as that pre- 
sented in Fig. 4. Moreover, it can be used even for the 
ground state (lowest energy for given j) giving a plausible 
fitting for very small and very large a, although serious 
discrepancies occur for a ~ 1 rendering it absolutely un- 
satisfactory in this region. 



(x 2n +i - x 2n ) 1 ~ 1 - S 2n + 5l n - <5f„. From Eqs.([35|) 
and (|36[) we get the nonlinear set of equations 



d 2 S 2n 
dr 2 



vf(S2n+l + fan-l - 2<5 2 „) 
(S 2n +1 + fon-l) ~ 4<5 2 n-l 



-4(n + l)(<5 2 2 „+i " - MSL-i ~ SL-i) > (37) 



C. Method of level dynamics 

The formalism of the generalized Calogero-Moser gas 
has been developed 0, 0, H3] as a useful alternative 
approach for Hamiltonians of the type H = Hq + aV 
[the perturbation term linearly dependent on a; Eq.© 
for a = p is the case with the first two terms as Hq 
and the term with a as V] by mapping the eigenvalue 
equations on a statistical many-body system of interact- 
ing pseudoparticles. There a is a controlling parame- 
ter considered as a pseudotime, a = r. The eigenval- 
ues were defined as dynamic variables E n (a) = x n (r), 
dxn/dr = Pn(r). The respective classical dynamic equa- 
tions for the pseudoparticles evolving in the pseudotime 
t are 



dp n 
dr 



E 



{Xyyi X n ^ 



dL r , 



dr 



where 



E 



LjniLin 



XI 



(33) 



d 2 s 



2n+l 



+45 2n+2 - 4(5 2 n 



+2 



dr 2 

= -/ + 4 

)-8|j- 

b" - i/2l), t 



v 2 (S 2n +2 + S 2n - 25 2n+ i) 



1 

~ 2 



(S 2n + s 2n+2 ) 



where v\ = 4(n + 1 + \j - 1/2|), v\ = 4(n + 1), and 
/ = 4 + 8|j-l/2|. 

For large n > \j\ and vf » w|,Eqs. (J3TJ) and ([38]) 
exhibit wavelike character in n. The linear approxima- 
tion to (|3"5j) bears the solution #2n+i ~ Jo(2\/|2j — l|n) 
where Jq(z) is the Bessel function. This behavior can 
be recognized in the upper part of the spectra in Figs. 
1(a) and 4(a). It is equivalent to the weak-coupling limit 
of the preceding subsection. The level avoidings in this 
linear part are kinematical and lead to clustering (dimer- 
ization) visualized in the level-spacing probability distri- 
butions (Sec. IV). 

Equations (f3"T)) and (|3"5|) allow a useful insight into 
the effects of nonlinearities on the level dynamics in- 
volved. Exploiting their wavelike character we define the 
phase variable £ = n — vt (v being an arbitrary con- 
stant velocity) and write an approximate equation for 
S 2n +i = $2n+i - 1 at \j\ > n: 



Pn = (n(T)\V\n(r)) = V n 



i(t) = {x n {r) - x m (r)) ■ V„, 



(34) 



The system of differential equations in the form |33|) 
is equivalent to the initial set |TTJ) and (fT5|) and is repre- 
sented as a set of dynamical equation of motion for a gas 
subject to a repulsive 1/r 2 potential ("Calogero-Moser 
gas"). In our case from Eqs. (fT"5)) and (|16[) we get (we 
simplify the notation, n = n r ) 



L 2n2n +i — V 2 

-^2n-12n 



n+l + |j-l/2|-(x2n+l(T) 

-V2Vn+l ■ (x 2 „(r) - 



- x 2n (r)), 



The dynamic equations for a pair of even and odd coupled 
trajectories are then 



>2n2n+l 



■ L 



2n+12n 



dp 2ri 



dp 2 „ _ ^1 
dr {x 2n+ i - x 2n ) 3 

r ,L 2 n+12n+2 • ^2rt+22n+l 



+ 2 



i-i • L 2 



dT 



{x 2n +2 ~ X 2n+1 ) 3 



{X2n-1 - X 2n ) 
L 2 n+12n ' L 2n2n ^ 

(x 2n - x 2n+ i) 3 



J2n+1 



d( 2 



+ 16|j|(Wi+^„+i)«0. (39) 



If 4|j| < v 2 , then 



S 2 n+l{C - 



Co) = --cosh" 



(4|j|) 



1/2 



{v 2 _ m) l/2 



(C-Co) 



(40) 



There Co = 
translational 
dx 2n+ i/dn = 
solution 



n — vt is a constant restoring the 
invariance of the solution; identifying 
l+^2n+i we recover a soliton-shaped (kink) 



4|j|) 
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1/2 



■ tanh 



(4|j|) 



1/2 



r(C-Co) 



(u 2_ 4^1)1/2' 

C - Co = n - n - v(t - r )(41) 



For what 

by x 2n+ i 



follows 

- x 2n 



define 
1 + 



a 

S 2n 



small fluctuation 
and approximate 



3 corresponding to the tunneling between two adjacent lev- 
j els. There, a multitude of true dynamical level avoidings 
-(36) r esults in a system of kink trains between subsequent 
pairs of levels. The limiting velocity |u| = 2y/J corre- 
5 sponds to the limit of vanishing fluctuation ^n+i an( i 
thus to suppressing the tunneling. 
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To consider now a situation interpolating between the 
two previous situations we take #2n = &2n + a; in Eq. ([55]) 
for §2 n , therefore choosing x so as to eliminate the 
quadratic terms in Eq. (|37[) . Let us take approximately 
<$2n+i ~ $2n— l ~ — S 2n . Then we arrive at equation 



1 d 2 8 2n d 2 8 2n 



4 dr 2 



dn 2 



In- 2\j\ ) S 2n +2nP 2n = i(2\j\-n/9) 



(42) 

If the oscillations in n can be neglected, under condi- 
tions for the coefficients of Eq.(j42]), B = 2n/3 - 2\j\ > 
and F ee (4/3)(2|j| - n/9) > (or n > 3\j\ > n/6), 
there exists an exact solution [2(| to (|l2"j) of periodic non- 
sinusoidal form spanned by F which describes the kink 
nucleation fluctuation. Increasing j causes the growth of 
the energy of the nucleus until a new kink generates. The 
respective solution reads 



S 2 u{t) 



n\ + cos(2u>r) 
n 2 + cos(2wt) ' 



(43) 



where a, w, n\, and n 2 are given as 



level. As already mentioned in the Introduction, the clas- 
sical counterpart of the system under consideration can- 
not be defined uniquely. The semiclassical approximation 
in two-level systems generically leads to classical chaotic 
patterns as a result of the nonlinear coupling between two 
subsystems, boson and electron, considered, respectively, 
as classical and quantum. The onset of classical chaos 
corresponds to energies above the first diabatic line. Re- 
ferring to Eq. (|46|) it is to be emphasized that the chaotic 
behavior refers to mentioned a purely quantum regime of 
medium values of j and n r between the weak coupling 
with j <C n r (dimcrized pairs of oscillators) and strong 
coupling with j 3> n r (kink lattice) domains. Thus, this 
chaotic behavior can be regarded as being of essentially 
quantum nature. One can conclude that the mapping of 
a quantum system onto the classical Calogero-Moser gas 
with repulsive interactions enables one to use the classi- 
cal formalism for describing the system via its quantum 
numbers. So a promising feature of this approach is its 
ability to represent quantum chaos by means of classical 
equations. 



ni = (2 - n\)/n 2 , 



B 



2B- 



2n2 + n\ 1 



2B 



(2 + ni). 



(44) 



From the relation F — a(2na 2 — B) it is evident that 
the amplitude a of the fluctuation S 2n is spanned by the 
driving field F, growing with \ j\ until the energy of the 
fluctuation reaches the energy of the kink. 

Equations (|37f and (|38|) also admit a chaotic regime 
for sufficiently large j,n under conditions d 2 8 2n /dT 2 = 
8\j - 1/2| + ci and d 2 S 2n+1 /dr 2 = -8\j - 1/2| + c 2 , 
where c\ = c 2 = 0, if the oscillations in n can be 
neglected- i.e. S 2n+1 + <5 2n _i - 25 2n w d 2 5 2n /dn 2 = 
and the same for 5 2n +i. One obtains then 5 2n + 5 2n + 2 = 
2n 2 8 2n+1 {l-8 2n+1 )/[\j~l/2\{n+j -1/2% where S 2n+1 = 
\j - l/2\5 2n+ i/n. Since 5 2n+2 « 5 2n , the famous logistic 
equation is recovered 



J2n+2 



AS 



2n+l 



(1 - $2n+l), 



(45) 



where A = n 2 /[\j - l/2|(n + j - 1/2)]. Equation (gSJ) 
assumes the transition to the chaotic regime for A > 
A crit = 3.56994. . ., i.e., for 



n C rit = \j—l/2\ 



Ac 



A r 



-- 2.1921-|2j-l|. 

(46) 

The condition n > n cr it specifies the values of n = n r 
from the nonlinear region above the kink nucleation 
threshold n > 3\j\ as estimated above. On the other 
side the maximum value of A, A = 4, implies, — 
(2j — 1)(1 + V2), which determines the upper limit for n, 
ii < 714. 

Traditionally the term "quantum chaos" is used to de- 
note the traces of classical chaotic behavior at a quantum 



IV. LEVEL STATISTICS 

Statistical methods of investigation of the complex 
spectra are based on evaluation of various "chaoticity 
degrees" of quantum systems - i.e., the density of levels, 
distribution of nearest-neighbor spacings (NNS's), and 
statistics of distribution of "curvatures" accounting for 
avoided crossings as well as more complicated measures 
of the distribution of avoided crossings 0, H3, HI] . Most 
of these characteristics can be obtained within the frame- 
work of random-matrix theory (RMT) [U [H \M, H3, 
whose predictions are supposed to adequately describe 
the quantum systems with generically chaotic classical 
analogs. However, the predictions of RMT cover in a 
satisfactory fashion only the " chaotic" limit of these sys- 
tems, classifying them into three universality classes ac- 
cording to the time-reversibility of the underlying Hamil- 
tonian [53, H^, [3l[ (that is, Gaussian orthogonal GOE, 
unitary GUE, and symplectic GSE ensembles). The 
fourth universality class is sometimes marked as corre- 
sponding to regular dynamics- for example, if the clas- 
sical analog of the system for some values of controlling 
parameters exhibits mostly regular KAM trajectories. 

In recent years, however, there came an understand- 
ing that chaotic patterns are to be sought not only from 
the energy-level distribution themselves but rather from 
pictures of their changing with the change of controlling 
parameters in the Hamiltonian (a). Thus, an alternative 
statistical approach is based on the concept of fully inte- 
grable description of the spectra by the dynamical model 
of a Calogero-Moser gas of interacting (repulsing) pseu- 
doparticles plj [2(| where the controlling parameter is 
considered a pseudo-time (Section IIIC). The dynamical 
equations for the eigenfunctions |s;(t)) of this Hamilto- 
nian bring, for example, to the notion of level curvature 
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FIG. 5: The widths of curvature distributions AK = (A 2 K) as function of j (a) and a (b). The scaling laws AK ~ j" and 
Aif ~ q m are seen. 



x(t) [13, HI] and several conjectures about its relation 
to the nature of the underlying RMT ensemble. Various 
characteristics of the level curvature distributions were 
reported as handy indicators of quantum chaos traces in 
a system. In particular, it was suggested [2l| that for 
energy regions with quantum chaos there holds the scal- 
ing AK ~ j" with positive v; meanwhile, for nonchaotic 
regions v is rather negative or close to zero. We inves- 
tigated the widths AK of the distributions of the level 
curvatures P(K) displayed as functions of j and a [Fig. 
5(a), 5(b)]. The level statistics of Fig. 5 (and of the fol- 
lowing Fig. 6 and 7) was taken from the energy intervals 
above the first diabatic line where multiple avoided cross- 
ings come into being. For the improvement of statistics 
we used the standard recipe of collecting the level se- 
quences from a number of a values lying in the neighbor- 
hood. Fig. 5(a) apparently shows the scaling behavior 
AK ~ j" with the scaling coefficient v G (0.99 - 1.01) 
for all a. In Fig. 5(b) we observe an another interest- 
ing regularity: namely, the scaling law AK ~ a M with 
/i ~ 0.5 ± 0.05. This latter observation opens a challeng- 
ing suggestion about the diffusive character of spreading 
the probability function P(K) in pseudotime a. From 
a closer inspection of these distribution functions one 
can even anticipate that the resulting diffusive equation 
would have a telegraph term and possibly is capable of 
bearing anomalous diffusion patterns. We are going to 
return to this point elsewhere. 

The realm of the symmetric JT model with definite j 
value is that of an effectively two-dimensional oscillator 



(or rather two oscillators pertaining to two electron lev- 
els); thus, the level-spacing picture for the unperturbed 
problem a — is trivial: all levels of a quantum harmonic 
oscillator are equidistant, so P(S) = S(S — 1). The non- 
trivial behavior of the NNS distribution appears to be 
consequence of the effective coupling of two oscillators 
through a. In Fig. 6 we show the level-spacing distri- 
butions P(S) for a set of j with different a. Our results 
indicate that the most close to RMT prediction situations 
[RMT would yield a Wigner distribution of NNS's in the 
form P(S) ~ Sexp(— S 2 )] are located at intermediate 
values of j and a as it is seen for pictures with j = 63/2. 
Meanwhile, the extreme values of j and a do not fit RMT 
at all. However, one can notice a crossover between two 
distinctly different types of behavior seen, for example, 
at j = 63/2, a = 2 and j = 255/2. Such a behavior 
can be semiquantitatively explained on the basis of the 
conception of effective potential wells- for example, built 
on the semiclassical coherent probe function in the space 
of their parameters (see 18]; this variant of building an 
effective potential is in fact the Husimi representation of 
the Hamiltonian operator taken for zero values of clas- 
sical momenta). Namely, the characteristic for big j is 
the presence of two almost noninteracting wells which are 
seen from the right-hand of Fig. 4 showing spectral en- 
tropies. We can assume that each of the wells refers to a 
separate quantum oscillator having level spacing a and b 
(let a < b). The superposition of the distribution for two 
independent sets of levels is found to be in the form (see 
for details, e.g., Berry et al(32j): 
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P(S) a=2, j=63/2 P(S) a=5, j=63/2 




FIG. 6: (Color online) Level spacing distributions lor particular j. Dimerization clustering related to two dominating level 
spacings is evident for a — 2 and j = 63/2. The third peak from localized states occurs at large j and a, dominating close to 
the semiclassical limit (compare with Fig. 4). The crossover from the dimerized regime to the localized regime with increasing 
j and a is evident. 



P{S)~~9(a-S)+5(S-a), (47) 
o 

9 being step function [we considered a <C b; if they 
are comparable, a refinement to (|47p is straightforward]. 
Conformity of Fig. 6 to the prediction (|47|) with a ~ 1.7 
and 6 = 4 for large j and a (for j = 255/2, a > 5) 
is perceivable. From the "adiabatic sheet" treatment in 
fl8| one can even define the parameters of both effective 
potential wells. The opposite case of small j (a) appar- 
ently indicates a dimerized structure with two dominat- 
ing spacings, with increasing extent of correlation when 
decreasing j and a. The graphs of the M-shaped proba- 
bility distributions for fixed j in Fig. 6(a) for small j and 
large n are similar to those by Cibils et al. [H| for s = 1/2 
which corresponds to our parity p = 1. There, the asym- 
metry between the matrix elements f nn and f n +i,n is 
negligible and the weak-coupling approach is satisfactory. 
However, for large j with large asymmetry of the matrix 
elements the strong-coupling limit of kink lattice is valid. 
There the relief of the three-minimum ground-state po- 
tential unfolds (Figs. 4) and the three-peak picture arises 
[Figs. 6(b) and 6(c)]. 

The level-spacing statistics taken separately for par- 
ticular values of j is the basis of a more realistic inves- 
tigation of the level statistics of the whole spectrum of 
the system whose Hamiltonian in this representation is 
a block matrix of blocks on the diagonal corresponding 
to definite j and off-diagonal blocks coupling [in the gen- 
eralized Ecg)(&i + b 2 ) JT model] for different j's. The 
peculiarity of the symmetric Jahn- Teller system is that 
the interblock coupling is exactly zero, and as was shown 
already in [33i ] with a variant of the central limit theo- 
rem, the superposition of (infinite or very large number) 
independent blocks will give the cumulative distribution 
which must be of Poissonian character, P(S) = exp(— S). 
This result is exactly that obtained by Yamasaki et al. 
[T3] for the absence of (trigonal) nonlinearity and which 
was claimed to be testimony of the absence of any traces 



of quantum chaos in this case. The detailed account for 
the level statistics for the generalized JT model includ- 
ing the symmetric JT as its partial case is reported by us 
elsewhere [HJ • For the sake of the present paper we just 
mention the serious deviations of the expected Poisson 
result in the symmetric case in the domain of dominat- 
ing quantum fluctuations (1 > (a,/?)). Indeed, the P(S) 
distribution shown for a — (3 — 0.5 exhibits consider- 
able deviations from the expected Poisson shape, show- 
ing anomalously large variance > 4 (for a = f3 — 1 this 
dispersion is ~ 1.17, and for the domain 1 < a it is rather 
rigid, falling normally within 0.7 — 0.8 irrespectively to 
the values of a). These anomalies are apparently due to 
the high nonuniformity of the spectral statistics in this 
case for different energy intervals. In Fig. 7 we show the 
level-spacing distributions for a — [3 — 0.5 taken for sep- 
arate level intervals (100 - 500, 500 - 1000, 1000 - 1500, 
and 1500— 2000) from which this difference is apparent. 
The distributions with a, (3 > 1, in contrast, are rather 
robust with respect to changing energy intervals. 



V. CONCLUSION 

The energy spectrum of a long but finite chain of corre- 
lated clusters is governed by the interplay (competition) 
of intracluster and intercluster interactions (repulsions) 
determined by two quantum numbers: rotational j and 
radial n r . At small \j\ <C n r , the excited spectrum 
is long-range-ordered and the broad dimerized region is 
a sample of pairwise oscillating levels. Respective level 
avoidings in the upper part of the spectra [the inset in 
Fig. 1(a)] result from the oscillations due to comparable 
intracluster and intercluster interactions in this part of 
the spectra. Such a regular wavelike spectrum is known 
to be characteristic for one-boson two-level systems since 
long time ago [35| and is a signature of the dimerization. 
The similarity of this part of the E®e spectra with two 
bosons to the spectra of one-boson systems was discussed 
in a more detail in the Introduction. 
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In the strong coupling limit of large j, \j\ 3> n r (the 
intercluster interaction is small or negligible when com- 
pared to the intracluster one) there appears a kink lat- 
tice - i.e., the regular lattice of oscillating level clusters. 
There the levels are bridged by a flip (kink) up to a higher 
level due to the tunneling [Fig. 1(b) above the diabatic 
line and inset therein]. At moderate values of \j\ and a 
there appears an intermediate region of strong nonlinear 
fluctuations responsible for the nucleation of the kinks 
[Figs. 1, 4(c),and 4(d)]. 

By an approximate quantum treatment of both the 
strong- and weak-coupling limit in the Sec. Ill we have 
found analytically the respective wavefunctions and en- 
ergies of localized exotic and extended states in a form 
of modified coherent oscillators: at strong coupling we 
showed the appearance of two branches of coherent os- 
cillators with effective frequencies f^// = 2 ± a/ \/2j — 1 
squeezed or antisqueezed by j to the degenerate quan- 
tum limit = 2 at j — > oo. An excellent agreement of 
these solutions to the exact numerical wave functions for 
large j (> 63/2) was marked. In the weak-coupling limit 
of small j <C n r the first order of perturbation theory 
yields satisfactory agreement with exact results for the 
wave functions except for a ~ 1 (f2 ~ 1). Therefore sim- 
ple approximate analytical solutions for wave functions 
are available in both the strong- and weak-coupling limits 
corresponding to the soliton (kink) and wavelike lattice 
(in the plane energy-coupling parameter). 

The mapping of our model on the Calogero-Moser gas 
of classical pseudoparticles with repulsive interactions en- 
abled us to describe (Sec. Ill C) all regions of the spectra 
by a set of classical equations in terms of its quantum 
numbers. In the intermediate region the equations under 
certain conditions implied the logistic equation route to 
chaos in terms of quantum numbers. 

On the other hand, the qualitative features of the ex- 
cited spectra of the E(g>e Jahn- Teller model can be un- 
derstood from the outline of the shape of the effective 
potential used for investigating the ground state of the 
model [H, H3] • With increasing j and a up to the semi- 
classical limit there emerge up to three effective potential 
minima. In accordance with this we identify the above- 
mentioned three phases in the spectra: for small coupling 
a the dimerized phase is related to two broad minima of 
the potential. Increasing a opens one additional nar- 
row minimum (for the ground state its counterpart was 
responsible for emerging the "tunneling phase" or light 
polaron [36|). The complex interplay of wave functions 
located over three minima corresponds to our intermedi- 
ate region [Fig. 4(c) and 4(d)] where the level statistics 
(Fig. 6) exhibits chaotic patterns close below the limit 
expected for fully chaotic RMT statistics [29(. Next, sub- 
sequently increasing a and j brings a suppression of the 
initial two "weak-coupling" wide minima in the quan- 
tum limit of large j. Thus there persist two almost non- 
interacting potential wells: a narrow and a wide one. 
Two remaining noninteracting branches of the spectra 
are markedly seen in Fig. 4(f). The lowest branch of the 



spectral entropy in Fig. 4(f) is related to the (third) nar- 
row minimum of the effective potential responsible for the 
kink lattice phase of avoided crossings (localized states) 
in energy spectra. This behavior is also supported by 
the shape of the level-spacing distribution [Fig. 6(d) and 
Eq. (|47p ] which has the marked character of superimpos- 
ing two almost noninteracting effective oscillators with 
different frequencies. 

The statistical level-spacing probability distributions 
investigated in Sec. IV were numerically calculated for 
particular j's. Their features are in agreement with above 
conclusions based on the shape of the effective potential: 
They are shown to exhibit up to three maxima (Fig. 6) 
depending on j and a. For moderate j and a a pair of 
almost symmetric prominent peaks of P(s) are related to 
two dominating strongly correlated potential wells which 
are getting closer at small a and j. The two dominating 
peaks of level spacings refer to the dimcrization cluster- 
ing of subsequent even and odd levels. At large j and 
a the NNS distributions are strictly determined by the 
shape of the ground-state effective potential (Fig. 6(d). 
The analysis of the cumulative statistical distributions of 
level spacings shows a strongly changing picture of differ- 
ent segments of spectra especially at weak couplings (Fig. 
7) . The dimerized phase is pronounced in the lowest part 
of spectra where a strong clustering is apparent from two 
peaks of the level-spacing distribution P(S). In the mid- 
dle part of spectra some patterns of the Wigner distri- 
bution at least at large S are perceivable. The widths 
of the level curvature distributions were found (Fig. 5) to 
be scaled as AK oc j", v oc 1, and thus satisfying the 
criterion suggested [2 11 ] for indicating chaotic patterns in 
our quantum system. Additionally, the scaling behavior 
AK oc a 5 was revealed which traces a possible bridge 
between mechanical and statistical points of understand- 
ing the complexity in the system referring to a possibility 
of reformulating the level curvature statistics stochasti- 
cally, as a Calogero-Moser gas diffusing with "time" a 

Since the sets of levels with j =/= j for the symmetric 
JT system are uncorrelated (matrix elements between the 
states are zero), the respective levels intersect (Sec. II). 
The cumulative level-spacing probability distributions 
consist of independent contributions of all j s and, con- 
sequently, the Poissonian distribution should supercede 
the above-described partial distributions. The cumula- 
tive level-spacing distributions, however, are found to ex- 
hibit enormously large dispersions 3> 1 at small a < 1 
(Fig. 7) which can be evidence for the overlapping of 
distributions of two kinds: a Poissonian one from the 
crossing of uncorrelated levels j ^ j and two-peak ones 
for the particular j's described above. In the quantum 
limit j — ► oo the spectrum tends to a continuum and the 
respective levels at j — > j will tend to avoiding levels. 
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FIG. 7: (Color online) Level-spacing distributions for subsequent segments of cumulative energy spectrum for small a < 1. 
Dimerization clustering peaks as remnants of the distribution in Fig. 6 are evident. The reference curves pertain to Wigner 
(solid line), semi-Poisson (dotted line), and Poisson (dashed line) distributions. 



